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The last couple of years has seen a remarkable number of new, explicit examples of the Wiener- 
Hopf factorization for Levy processes where previously there had been very few. We mention 
I— I in particular the many cases of spectrally negative Levy processes in [21, 32], hyper-exponential 

and generalized hyper-exponential Levy processes [24], Lamperti-stable processes in [9, 10, 13, 39], 
^ Hypergeometric processes in [35, 31, 11], /3-processes in [29] and 6'-processes in [30]. 

In this paper we introduce a new family of Levy processes, which we call Meromorphic Levy 
processes, or just M-processes for short, which overlaps with many of the aforementioned classes. 
C A key feature of the M-class is the identification of their Wiener-Hopf factors as rational functions 

of infinite degree written in terms of poles and roots of the Laplace exponent, all of which are real 
CN numbers. The specific structure of the M-class Wiener-Hopf factorization enables us to explicitly 

^ handle a comprehensive suite of fluctuation identities that concern first passage problems for finite 

and infinite intervals for both the process itself as well as the resulting process when it is reflected in 
its infimum. Such identities are of fundamental interest given their repeated occurrence in various 
fields of applied probability such as mathematical finance, insurance risk theory and queuing theory. 
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H 1 Introduction 



The theory of Levy processes forms the cornerstone of an enormous volume of mathematical literature 
which supports a wide variety of applied and theoretical stochastic models. One of the most obvious 
and fundamental problems that can be stated for a Levy process, particularly in relation to its role as 
a modelling tool, is the distributional characterization of the time at which a Levy process first exists 
either an infinite or finite interval together with its overshoot beyond the boundary of the interval. As 
a family of stochastic processes, Levy processes are now well understood and the exit problem has seen 
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many different approaches dating back to the 1960s. See for example [2, 6, 8, 16, 41, 19, 13, 35] to name 
but a few. 

Despite the maturity of this field of study it is surprising to note that, until very recently, there were 
less than a handful of examples for which explicit analytical detail concerning the first exit problem 
could be explored. Given the closeness in mathematical proximity of the first exit problem to the 
characterization of the Wiener-Hopf factorization, one might argue that the lack of concrete examples of 
the former was a consequence of the same being true for the latter. The landscape for both the Wiener- 
Hopf factorization problem and the first exit problem has changed quite rapidly in the last couple 
of years however with the discovery of a number of new, mathematically tractable families of Levy 
processes. We mention in particular the many cases of spectrally negative Levy processes in [21, 32], 
hyper-exponential and generalized hyper-exponential Levy processes [24], Lamperti-stable processes in 
[9, 10, 13], Hypergeometric processes in [35, 31], /3-processes in [29] and 6'-processes in [30]. 

In this paper we introduce a new family of Levy processes, which we call Meromorphic Levy processes, 
or just M-class for short, that overlaps with many of the aforementioned classes. Our definition of the 
M-class of processes will allow us to drive features of their Wiener-Hopf factors through to many of 
the fluctuation identities which are of pertinence for a wide variety of applications. In the theory of 
actuarial mathematics, the problem of first exit from a half-line is of fundamental interest with regard to 
the classical ruin problem and is typically studied within the context of an expected discounted penalty 
function. The latter is also known as the Gerber-Shiu function following the first article [18] of a long 
series of papers found within the actuarial literature. In the setting of financial mathematics, the first 
exit of a Levy process, as well as a Levy process reflected in its infimum, from an interval is of interest 
in the pricing of barrier options and American-type options, [1], as well as certain credit risk models, 
[20, 33]. In queueing theory exit problems for Levy processes play a central role in understanding the 
trajectory of the workload during busy periods as well as in relation to buffers, [17, 27]. Many optimal 
stopping strategies also turn out to boil down to first passage problems for both Levy processes and 
Levy processes reflected in their infimum; classic examples of which include McKean's optimal stopping 
problem, [38], as well as the Shepp-Shiryaev optimal stopping problem, [2, 3] . 

It is not our purpose however to dwell on these applications. As alluded to above, the main objective 
will be expose a comprehensive suite of fluctuation identities in explicit form for the M-class of Levy 
processes. We thus conclude the introduction with an overview of the paper. 

In the next section we give a formal definition of meromorphic Levy processes, dwelling in particular 
on their relationship with discrete completely monotone functions. In Section 3 we consider several 
classes of existing families of Levy processes that have appeared in recent literature. Next, in Section 
4 we establish explicit identities for the exponentially discounted first passage problem. In particular 
we deal with (what is known in the actuarial literature as) the Gerber-Shiu measure describing the 
discounted joint triple law of the overshoot, undershoot and undershoot of the maximum at first passage 
over a level as well as the marginal thereof which specifies the law of the discounted overshoot. It is 
important to note that the discounting factor which appears in all of the identities means that one is 
never more than a Fourier transform away from the naturally associated space-time identity in which 
the additional law of the time to first passage is specified. This last Fourier inversion appears to be 
virtually impossible to produce analytically within the current context, but the expressions we offer are 
not difficult to work with in conjunction with straightforward Fourier inversion algorithms. 

In Section 5 we look at the more complicated two-sided exit problem. Inspired by a technique of 
Rogozin [41] we solve a system of equations which characterize the discounted overshoot distribution 
on either side of the interval in question. The same technique also delivers explicit expressions for the 
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discounted entrance law into an interval. In Section 6 we look in analytical detail at what can be said of 
the ascending and descending ladder variables. In particular we offer expressions for their joint Laplace 
exponent and jump measure. Section 7 mentions some additional examples of fluctuation identities which 
enjoy explicit detail and finally in Section 8 we describe some numerical experiments in order to exhibit 
the genuine practical applicability of our method. 

2 Meromorphic Levy processes 

Recall that a general one-dimensional Levy process is a stochastic process issued from the origin with 
stationary and independent increments and almost sure right continuous paths. We write X = {Xt : 
t > 0} for its trajectory and P for its law. As X is necessarily a strong Markov process, for each x G M, 
we will need the probability P^: to denote the law of X when issued from x with the understanding 
that Pq = P. The law P of Levy processes is characterized by its one-time transition probabilities. In 
particular there always exist a triple (a, a, 11) where a G M, a G M and 11 is a measure on M\{0} satisfying 
the integrability condition J^{1 A x'^)Il{dx) < oo, such that, for all z G M 

E[e''^'] = e*^('^\ (1) 
where the Laplace exponent iplz) is given by the Levy-Khintchine formula 

^(z) = \a^z^ + fxz+ [ (e"^ - 1 - zxhix)) U{dx). (2) 

Here h{x) is the cutoff function which is usually taken to be h{x) = I{|a;|<i}. However, everywhere in the 
present paper we will work with the Levy measures H(dx) which have exponentially decaying tails, thus 
we will take h{x) = 1. Note, that the exponential decay of the tails of the Levy measure also implies that 
the Laplace exponent can be analytically continued into some vertical strip a < Ile{z) < b for a < < 6. 

Everywhere in this paper we will denote the tails of the Levy measure as Il~^{x) = H((a;, oo)) and 
n~(x) = H((— oo, — x)) for a; > 0. Let us define the supremum/infimum processes Xt = sup^^^Xs, 
2£_f = infsKtXs and denote by e(g) an independent, exponentially distributed random variable with rate 
g > 0. Finally, the first passage time above x is defined as = mf{t > : > x}, and similarly 
r- = inf{t >0:Xt<y}. 

Let us recall some basic facts about completely monotone functions. A function / : (0, oo) i— M is 
called completely monotone if / G C°° and (— l)"/*^"'^(x) > for all n G N U {0} and x > 0. According 
to Bernstein theorem (see Theorem 1.4 in [42]), the function / is completely monotone if and only if it 
can be represented as the Laplace transform of a positive measure on [0, oo): 

fix) = j e-"Xdz), X > 0. (3) 

[0,oo) 

Note that if /(O"'") = 1, then /z(dx) is a probability measure and 1 — /(x) is the cumulative distribution 
function of a positive infinitely divisible random variable Z, whose distribution is a mixture of exponential 
distributions. We will denote the class of completely monotone functions as CM.. 

Next, let us introduce a subclass of completely monotone functions, which will be important for us 
later. We will call / : (0, oo) i— > M a discrete completely monotone function if the measure /i(dz) in the 
representation (3) is discrete, and the support of the measure ^{dz) is infinite and does not have finite 
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accumulation points. This implies that the measure /i(dz) is an infinite mixture of atoms of size a„ at 
the points bn 

n>l 

where Sf,{dx) denotes the Dirac measure at x = b, for all n > 1 we have ctn > 0, &„ > and 6„ — )■ oo as 
n — 7- oo (without loss of generahty we can also assume that the sequence {&n}n>i is strictly increasing). 
From (3) it follows then that any discrete completely monotone function can be represented as an infinite 
series of exponential functions 

/(x) = J]a„e-''"^ x>0. (4) 

n>l 

We will denote the class of discrete completely monotone functions as VCAi. 

Definition 1 (M-class). A Levy process X is said to belong to the Meromorophic class (M-class) if 
n+(a;),n-(x) G VCM. 

We see that according to our definition of the discrete completely monotone functions, the process is 
Meromorphic if and only if the Levy measure n(da;) has a density with respect to the Lebesgue measure, 
given by 

7r(x) = I{^.>o} a„p„e-''"^ + I{^<o} ^ a^p^e'^"^, (5) 

n>l n>l 

where all the coefficients a„, a^, p„, p„ are positive, the sequences {pn}n>i and {p„}„>i are stricly 
increasing, and p„ — t- +oo and p„— T-H-ooasn— t-H-oo. 

Proposition 1. Assume that tt{x) is given by (5). The integral JjgX^7r(x)dx converges if and only if 
both series X]n>i ^nPn^ and X]n>i ^nPn"^ converge. 

Proof. The "if part was established in Proposition 1 in [30]. The "only if part follows from (5) and 
the Monotone Convergence Theorem. □ 

Recall, that a function g : C i— > CU{oo} is called meromorphic if it does not have any other 
singularities in the open complex plane except for poles. A function g{z) is called real meromorphic 
function if it is meromorphic and g{z) G M U {oo} for all z G M, or equivalently, if g{z) = g{z). 

Theorem 1. The following conditions are equivalent: 

(i) X is Meromorphic. 

(a) Tl~^{x),Tl~{x) G CM. and the Laplace exponent tp{z) is meromorphic. 

(Hi) For some q > (and then, for every q > 0) the functions F{Xe{q) > x) andF{—2Le(q) > ^) restricted 
to X > belong to the class VCM.. 

(iv) For some q > (and then, for every q > 0) the functions P(Xe(q) > x) andF{—Xc(g) > x) restricted 
to X > belong to the class VCM . 
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(v) For some q > (and then, for every q > 0) we have the factorization 



1 - ^ 1 + ^ 

g-m = <iUT^llY^^ zGC, (6) 

n>l Pn n>l Pn. 

where all roots {(n, —Cn} of iIj{z) — q are real and interlace with the poles {p„, —pn} 

• •• - P2< -6 < -pi < -Ci < < Ci < Pi < C2 < P2 < ••• (7) 

(vi) The Laplace exponent ipi^) is a real meromorphic function, which satisfies lm{ip{z)/z) > for all 
z in the half-plane lm{z) > 0. 

Proof. The main ideas and tools needed for the proof of this theorem come from the proofs of Theorem 
2 in [40] and Theorem 1 in [30]. 

(i) ^(ii) Using (5) and (2) we find that 

i^{z) = \a'z^ + pz + z'Y. , +^'Y. ^i^\v (8) 

2 ^ Pn[pn -Z) ^ Pn{pn + z) 

which shows that if'iz) is a meromorphic function. Since VCM. C CM. this proves (ii). 

(ii) ^(vi) We know that ip{z) is a meromorphic function, and that ip{'S) = 0, therefore ipiz) is analytic 
in some neighbourhood of zero. From the proof of Theorem 2 in [40] we find that (a) ip{z) can be 
analytically continued in the half-planes lm{z) > and Im(2;) < 0, (b) lm(il){z)/z) > for all z in the 
half-plane lm{z) > 0. Since ip{z) = ipiz) we conclude that ipiz) is a real meromorphic function. Using 
this fact and the above statement (b) we obtain (vi). 

(vi)^(v), Wq > 0. For g > and Im(2;) > it is true that lm{—q/z) > 0, thus (vi) implies that 
lm{{ip{z) — q)/z) > for all z in the half-plane lm{z) > 0. Using Theorem 1 on page 220 in [36] (the 
statement of this result can also be found in the proof of Theorem 1 in [30]) we find that (v) is valid for 
all g > 0. 

(v)^(iii) This follows from Theorem 1 in [30]. 

(iii) ^(iv) One can also check that if two functions /(x) and g{—x) belong to the class PCA^, then the 
same is true for I{a.>o}(/ * 9){.x) and I{x<o}(/ * 9){.^)^ where f * g is the convolution. The result (iv) then 

follows easily from this fact and the Wiener-Hopf decomposition which says that Xe(q) = Y1 + Y2, where 

the random variables Yi and Y2 are independent, Yi = Xe{q) and Y2 = 2Le{q)- 

(iv) ^(i) From (iv) and the definition of the T>CM. class we know that for some g > there exist cto > 
and positive constants d„, /3„ such that 



P(Xe(g) G dx) = ao5o(da;) + 



l{x>0} an^n + l{x<0} 



n>l n>l 



dx. (9) 



Note that the condition P(Xe(q) G M) = 1 implies that both series J2n>i ^nf^n ^ J2n>i ^n^n ^ converge. 
Let us define iIj{z) = z'^'E[exp{zXc{q))]. Using (9) we obtain 



4^{z) = aoz' + z'y ^ ^ +z'y.-^ , zee. 
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Comparing the above formula with (8) we conclude that iplz) is a Laplace exponent of a Meromorphic 
Levy process. We have already proved that (i) implies (vi), therefore lm{ilj{z)/z) > for all z in the 
half-plane lm{z) > 0. 

Next, using the definition of the Laplace exponent (1) it is easy to verify that E[exp(zXe(5))] = 
q/{q ~ i'iz)). Therefore 

ip{z) — q qz 

z 4,{zy 

As we have already established, lm{ilj{z)/z) > for all z in the half-plane lm{z) > 0. Using this fact and 
the above identity we conclude that lm{{ilj{z) — q)/z) > 0. Applying Theorem 1 on page 197 in [14] (the 
statement of this result can also be found in the proof of Theorem 1 in [30]) we find that ip{z) admits a 
representation of the form (8), which in turn implies that the process X is meromorphic. 

Finally, note that if (v) is true for some q > 0, then (iii) and (iv) are valid for the same value of g > 0. 
But as we have already demonstrated, any of the conditions (iii), (iv), (v) implies (i) which is equivalent 
to (v) being valid for all g > 0. Thus, if one of the conditions (iii), (iv), (v) is valid for some q > 0, then 
it must be valid for all g > 0. □ 



Statement (ii) in Theorem 1 shows that the M-class of Levy processes might also be called "processes 
with completely-monotone Levy measure and meromorphic Laplace exponent", which explains the ori- 
gin of the name "Meromorphic Levy processes". Note, however, that there exist Levy processes with 
meromorphic Laplace exponent but not completely monotone Levy measure. 

The following technical result on partial fraction decomposition of infinite products will be very 
important for us later. 

Lemma 1. Assume that we have two increasing sequences p = {pn}n>i o-nd ( = {Cn}n>i of positive 
numbers, such that pn — )■ +oo as n +oo and the following interlacing condition is satisfied: 

(10) 



Define 



Then for all z > 



Ci < Pi < C2 < P2 < 
1 + 



Hi 



Pn 



2 ' 



z>0. 



n>l Cn 



ao(p,C) + X]a„(p,C) 



Cn 



n>l 



Cn + Z 



i + zho{(,p) + Y,MC,p) 



Pn + Z 



where 



ao(p,C) 



bo(C,p) 



Ck 

lim TT— , a„(p,C) 
'^+°° Pk 



Cn 
Pn 



Cn 
Pk 



k>l 
ky^n 



Cn 
Cfe 



Ci " 



lim n/^' b„(C,p) = - 1 



Hi 

k>l 
k^n 



1 - 



Pn 

Pk 



(11) 

(12) 
(13) 

(14) 
(15) 
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Moreover, ao(p, C) > 0, bo(C, p) > and for all n > 1 we have a„(p, C) > 0, b„(C, p) > 0. 

Proof. The convergence of the infinite product in (11) and the partial fraction decomposition (12) follow 
from Theorem 1 in [30]. To prove the second partial fraction decomposition, rewrite the infinite product 
as 




and note that sequences {Cn+i}n>i and {p„}„>i satisfy interlacing condition, thus we can apply the same 
method. The details are left to the reader. □ 



Using Monotone Convergence Theorem one can show that formulas (12) and (13) are equivalent to 



ao(p,C) + 



X^a„(p,C)Cne" 



n>l 



-""dx. 



^b„(C, p)Pne ''"^ 



n>l 



(1 



e ^'') dx. 



(16) 
(17) 



Before stating our next result we recall that a Levy process creeps upwards if for some (and then 
all) X > 0, ^{Xr+ — ^\^x < ^) > 0- Moreover, we say that is irregular for (0, cxd) if and only if 
F{tq > 0) = 1. (Note that this probability can only be or 1 thanks to the Blumenthal zero-one law). 
We refer to Bertoin [5], Doney [15] or Kyprianou [34] for more extensive discussion of these subtle path 
properties. 

We also need to introduce some more notation. In the forthcoming text everything will depend on 
the coefficients {a„(p, C), a„(p, C)}n>o defined using (14) and {bn(C, p), b„(C, p)}„>o defined using (15). 
We define for convenience a column vector 

a(p, C) = [ao(p, C), ai(p, C), a2(p, C), -f , 

and similarly for a(p, Q, b(C, p) and b(C, p). Next, given a sequence of positive numbers ( = {Cn}n>i, we 
define a column vector y{(, vector of distributions 

v(C,x)= [5o(a:),Cle-«^^C2e-^^^...]^, 
where ^0(2;) is the Dirac delta function at x = 0. 

Prom here on, unless otherwise stated, we shall always assume that X is a Levy process 
belonging to the M-cIass and that X is not a compound Poisson process. 

Theorem 2 (Properties of Meromorphic processes). 



(i) The Wiener- Hopf factors are given by 

1 + 



iz) = E 



nr 



Pn 



- + -' 

n>l Cn 



E [e^^^(«)] = n I 



1 + 



n>l 



+ 
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(ii) For X > 

P(Xe(,) e dx) = a(p, C)'^ X v(C, x)dx, P(-X,(^) e dx) = a(p, C)'^ X v(C, x)dx. (19) 

(in) ao(p,C) (ao(p, 0) nonzero if and only ifO is irregular for (0, oo) (correspondingly, (— oo,0)). 

(iv) bo(CjP) (bo(C)P)) nonzero if and only if the process X creeps upwards (correspondingly, down- 
wards) . 



(v) For every q > 



P(Xe(g) edx) = q 



n> 



dx. (20) 



Proof. Parts (i) and (ii) were established in Theorem 1 in [30], these results also follow easily from 
formulas (6), (16) and the structure of the Wiener-Hopf factorization. Let us prove (iii). We note that 
is irregular for (0, oo) if and only if, for any q > 0, Xc(q) has an atom in its distribution at 0. From part 
(ii) this is clearly the case if and only if ao(p, C) is non-zero. For part (iv) of the proof, we note that we 
may necessarily write 

where K{q, z) is the Laplace exponent of the bivariate subordinator which describes the ascending ladder 
process of X. It is also known that a Levy process creeps upwards if and only if the subordinator 
describing its ladder height process has a linear drift component. The drift coefficient is then described 
by lirsiz^foo z) / z G [0,oo) where the limit is independent of the value of g > 0. Inspecting (17) we 
see that the required drift coefficient is non-zero if and only if bo(C5P) is non-zero. The conclusion thus 
follows. 

Finally, let us prove (v). From part (iv) of Theorem 1 we know that P(Xe(g) G dx) can be written in 
the form (9), where = Q due to our assumption that X is not a compound Poisson process. Combining 
this fact and the identity E[exp(2;Xo(g))] = q/{q — 'ip{z)) we conclude that 

= MPn- Z)^Y ^) ' -^C. (21) 

We see that the function on the right-hand side of the above equation has poles at the points /3„ and 
—Pn, while the left hand side has poles at (n and —(n (this follows from (6)). Therefore we conclude 
that (3n = Cn and /3„ = Cn- Comparing the residues of both sides of (21) at the pole z = we see that 

-Res ( ^ : z = (r, 



Similarly we find 



Res ( ^— - : z = -Cn ^ 



Combining the above two identities and (9) we obtain (20). □ 
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We would like to emphasize the importance of statement (v) of Theorem 2. It is well known that 
there exist only a few specific examples of Levy processes (such as Variance Gamma or Normal Inverse 
Gaussian processes) for which the law of Xf is known explicitly (for every t > 0). While we do not know 
the law of Xt for Meromorphic processes, the result in Theorem 2 (v) shows that Meromorphic processes 
have an advantage that at least the distribution of Xe(g) can be easily computed (for every g > 0). The 
formula is given in terms of the roots (n and (n, but as we will see in Section 8, computing these numbers 
is a rather simple task and it can be done very efficiently. 

As the next corollary shows, the Levy measure of a Meromorphic process X can be easily reconstructed 
from zeros and poles of q — ip{z). This has the spirit of Vigon's theory of philanthropy (cf. [43]) in that 
we construct the Levy measure from the Wiener- Hopf factors. 

Corollary 1. Assume thatq > and the Wiener-Hopf factors are given by (18). Then X is Meromorphic 
and it's Levy measure is given by (5), where {pn\n>i and {pn}n>i are poles of the Wiener-Hopf factors 
0g(— iz) and (pqiiz) and coefficients an and a„ are given by 

an = K{C,p)—^. — r, a„ = b„(C,p)— . (22) 

Proof. The fact that X is Meromorphic was already established in part (iii) of Theorem 1. Assume that 
the density of the Levy meausre 7r{x) is given by (5), then as we have established in the proof of Theorem 
1, the Laplace exponent ipi^) can be expressed in the form (8). Comparing this equation with (6) we 
conclude that the exponents p„ and p„ in (5) must coincide with the poles of ipi^z). 

From (8) it also follows that a„p„ = —Iies{ip{z) : z = Pn)- Using the Wiener-Hopf factorization 



q-ij{z)= " . (23) 

0g(-i2:)05(-i2) 

and the fact that (f)~{—iz) is non-zero at p„ we see that 

ttnpn = Res(0+(-i2;)"^ : z = Pn) X , , . — r- 

In order to finish the proof we only need to check that Res(0+(— iz)~^ : z = pn) = Pnb„(C,p) which 
follows from (13). □ 



Note that the coefficients and a„ given in (22) (which define the Levy measure via (5)) depend on 
q on the right hand side but not on the left hand side. It is therefore tempting to take limits as g J, 
on the right hand side. This is not as straightforward as it seems as the limits of a„ and d„ are not easy 
to compute. It would be more straightforward to start with the Wiener-Hopf factorization for the case 
that g = and then perform the same analysis as in Corollary 1. The next corollary provides us with 
the aforementioned Wiener-Hopf factorization. We recall that ^„ = Cn{,q) (resp. C„ = Cn{,q)) denote the 
non-negative solutions to equation ^/^(z) = q (resp. ip^—z) = q). 

Corollary 2. 

(i) Assume that E[Xi] > 0. As q we have q/Ci{q) — t- ]E[Xi] and for n > 1 

Cn+l(g)^Cn+l(0)7^0, Cn(g) ^ Cn(0) 7^ 0. 
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We have the Wiener- Hopf factorization —iplz) = k{0, —z)k{0, z) , where 



1 ^ > n>l Pn 

1 + 



n>l Pn 



('wj Assume that E[Xi] = 0. g ^ 0+ w;e have ^/Ci{q) ^ v^Var(Xi)/2, v^/Ci(g) ^ v^Var(Xi)/2 
anc? /or n > 2 

Cn(g) ^ Cn(o) ^ 0, Cn(g) ^ Cn(0) ^ 0. 
We have the Wiener-Hopf factorization —ipiz) = k{0, —z)k{0, z) , where 



9 ^ ^ n>l 



Proof. Let us prove part (i). Since E[Xi] > we know that Xc(q) — ?■ +oo and Xj,^^-) — )■ as 
g — )■ 0~^, and also that X^^ < oo with probabihty one. Therefore, if z > the Wiener-Hopf factor 
4>g{iz) = E[exp(— zXo(g))] must converge to as g — 0"*" while (f)~{—iz) = K[exp{z2(_c{q))] must converge 
to E[exp(zX^)] as g — )■ 0+. Since all roots (nil) and (nil) have nonzero limits as g — )■ 0"*" if n > 2 (this is 
true due to the interlacing condition (7) and the fact that ip{z) is a meromorphic function) we conclude 
that Ci{q) must go to zero while Ci(Q') ~^ Ci(0) 0. The function tp{z) is analytic in the neighbourhood 
or 2 = 0, and it's derivative ip'{0) = E[Xi] ^ 0, thus using the Implicit Function Theorem we conclude 
that Ci{q) is an analytic function of g in some neighbourhood of g = 0, and taking derivative with respect 
to g of the equation ip{z) = q we find that 

:Ci(g) 



dg'^^'^ ^'(Ci(g))' 

Since ip'{0) = E[Xi] we conclude that g/Ci(?) E[Xi] as g — )■ 0"*". Using this fact and the formulae in 
(18) for the Wiener-Hopf factors we obtain expressions (24). 

The proof of part (ii) is very similar. We use the fact that ^"'(0) = E[Xi] = and ip"{0) = Var(Xi) 
to conclude that ipi^z) = |Var(Xi)z^ + 0{z^) as z — > 0. This implies that when g is small and positive, 
the equation ^Ij{z) = q has two solutions in the neighbourhood of zero, Ci = A/2g/Var(Xi) + o{y/q) and 
— (^1 = — A/2g/Var(Xi) + o{y/q). Using this information, the uniquenes of the Wiener-Hopf factorization 
and taking the limit as g — i- 0^ in (6) we obtain (25). □ 



3 Examples of Meromorphic processes 

There are four particular families of Levy processes which have featured in recent literature, all of which 
have proved to have relevance to a number of applied probability models, all of which have exhibited 
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some degree of mathematical tractability in the context that they have occurred. These processes are as 
follows. 

• Hyper-exponential Levy processes: The elementary but classical Kou model, [28], consists of a 
linear Brownian motion plus a compound Poisson process with two-sided exponentially distributed 
jumps. The natural generalization of this model (which therefore includes the Kou model) is 
the hyper-exponential Levy process for which the exponentially distributed jumps are replaced by 
hyper-exponentially distributed jumps. That is to say, the density of the Levy measure is written 

N N 

7i{x) = I{^>o} a^piC-P''' + I{^<o} ^ aiPief^'"", 

i=l i=l 

where Oj, di, pi and pi are positive numbers and A^, N are positive integers. One can verify that 
the Laplace exponent is a rational function of the form (8), where we have finite sums instead of 
infinite series, and that hyper-exponential Levy processes have finite activity jumps and paths of 
bounded variation unless o"^ > 0. 

Strictly speaking, the Hyper-exponential Levy processes do not belong to the class of Meromorphic 
processes, as in this case the Laplace exponent ipi^) has only finite number of poles pn, —pn- 
However, all the results presented in this paper are still correct when we have a process with 
positive and/or negative hyper-exponential jumps, provided that the results are interpreted in the 
correct way, see Remark 1 below. 

• /3-class Levy processes: The /3-class of Levy processes was introduced by Kuznetsov [29]. The 
Laplace exponent is given by 

^{z) = V^2 + a^+^|B(ai-^,l-Ai)-B(ai,l-Ai)| 

+ 1 |B(a2 + ^, 1 - A2) - B(a2, 1 - As)} , 

where B{x,y) = r(x)r(?/)/r(x + y) is the Beta function, with parameter range a G M, cr > 0, 
Ci,ai, (3i > and Aj G (0,3) \ {1,2}. The corresponding Levy measure H(dx) has density 

With the help of binomial series one can check that 7t{x) is of the form (5), thus /3-processes are 
Meromorphic. 

The large number of parameters allows one to choose Levy processes within the /3-class that have 
paths that are both of unbounded variation (when at least one of the conditions cr 7^ 0, Ai G (2, 3) 
or A2 G (2, 3) holds) and bounded variation (when all of the conditions a = 0, Ai G (0, 2) and 
A2 G (0,2) hold) as well as having infinite and finite activity in the jump component (accordingly 
as both Ai, A2 G (1,3) or Ai, A2 ^ (1,3)). The /3-class of Levy processes includes another recently 
introduced family of Levy processes known as Lamperti-stable processes, cf. [9, 10, 13]. 

• ^^-class Levy processes: 
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The 6'-class of Levy processes was also introduced by Kuznetsov [30] as a family of processes with 
Gaussian component and two-sided jumps, characterized by the density of the Levy measure 

7r(a;) = I{.>o}Ci/3ie-'^^"e,(a;/3i) + I{.<o}C2/32e"^"efc(-x/32), 

where Ofc(x) is the k-th order (fractional) derivative of the theta function ^3(0,6^^): 

6^(0;) = ^^3(0, = I|,=o} + 2 ^ n''e--'^, x > 0. 

n>l 

The parameter x = A; + 1/2 corresponds to the exponent of the singularity of tt{x) at x = 0, and 
when X ^ {1/2, 3/2, 5/2} {x G {1, 2}) the Laplace exponent of X is given in terms of trigonometric 
(digamma) functions. Again, due to the fact that the density of the Levy measure is of the form 
(5) we conclude that 6'-processes are Meromorphic. 

• Hypergeometric Levy processes: These processes were introduced in Kyprianou et al. [35, 31] 
as an example of how to use Vigon's theory of philanthropy. Their Laplace exponent is given by 

V'(z) = ^aV + az + $i(-2)$2(2), (26) 

where a, a G M and $1, $2 are the Laplace exponents of two subordinators from the /3-class of Levy 
processes. Note that such Laplace exponents necessarily take the form 

<l>ie) = k + (50 + -^{B (1 - « + 7, -7) - B (1 - « + 7 + 9/(3, -7) }. (27) 



Compared with the previous examples, it is less obvious that hypergeometric processes are Mero- 
morphic, however it is not hard to establish this result with the help of Theorem 1. We know that 
/3-processes are meromorphic, thus the process X with the Laplace exponent ip{z) = ^i{—z)^2{z) 
is Meromorphic (this is due to part (v) of Theorem 1). Finally, the hypergeometric process defined 
by the Laplace exponent via (26) and the process X have the same Levy measure, therefore X is 
also Meromorphic. 

Remark 1. Within the scope of Definition 1, when the process Xt has hyper-exponential positive jumps 
(a mixture of N exponentials), we have only a finite sequence {pn}n=i,2,..,N and either N (or + 1) 
negative roots (n- All the formulas presented in this paper are still valid in this case if we adopt notation 
pk = 00 for k > N and Cfc = 00 for k > N {or k > N + 1). For example, expression (18) for the 
Wiener-Hopf factor becomes a finite product 

n=l Cn 

depending on whether we have A^ or A^-|- 1 negative roots (n- The same remark holds true when we have 
hyper-exponential negative jumps. 

Remark 2. When the process X has hyper-exponential positive jumps, calculation of coefficients a(p, Q 
and h{p,() is much easier, since the corresponding expressions in (14) and (15) are just finite products. 
Moreover, calculation of the coefficients a(p, C) and h{p,() is also simplified considerably even if the 



N 

n 



1 + 



pn 



1 1 + f 



1 + 



z ' 



CiV+1 
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sequence {p„}n>i is infinite, since we can use Wiener-Hopf factorization (23) and tlie fact tliat — C„ and 
—pn are simple roots (poles) of function iIj{z) — q. For example we can compute coefficients a„(p, C) as 
follows: 



an(p,C) 



-Res 



-iz) : z 



<r 



Cn 

—Res 

Cn V(?- V'W)0+(-i2;) 



-Cn. 



Cn^'(-C 



nyV^g (iCn) 



and this last expression is more convenient than (14), since usually we know il)'{z) explicitly and 0g (2) 
is just a rational function. 



4 One-sided exit problem 

Doney and Kyprianou [16] have given a detailed characterization of the one-sided exit problem through 
the so-called quintuple law. The latter can be easily integrated out to give the following discounted triple 
law which is also known in the actuarial mathematics literature as the Gerber-Shiu measure (cf. [7]). 
The following result holds for general Levy processes. 

Lemma 2. Fix c > 0. Define 

r+ = inf{t > : > c}. 

For all q,y,z > and m G [0, c V 2;] we have 



E 



'I(X_+ -cedy; c - X.+_ e dz ; c - X_+_ E du) 



P(Xe(g) G c - dM)P(-X,(^) edz- u)U{dy + z). 



Taking account of formulas (19) this gives us the following immediate corollary for the M-class of 
Levy processes. 

Corollary 3. Fix c > 0. For all q,y,z > and m G [0, c V we have 



E 



'QTc 



I(X^+ - cedy ; c - X^+_ e dz ; c- X^+_ G du) 



[a(p, C)^ X v(C, c — m)] a(/), C)^ X v(C, z — u) Il{dy + z)dzdu 



Often, one is only interested in the discounted overshoot distribution. In principle this can be 
obtained from the above formula by integrating out the variables z and u. However, it turns out to be 
more straightforward to prove this directly, particularly as otherwise it would require us to specify in 
more detail the Levy measure in terms of poles and roots for the M-class (cf. Corollary 1). The following 
result does precisely this, but in addition, it also gives us the probability of creeping. 



Theorem 3. Define a matrix A = {O'i,j}i,j>o as 

f 



a. 



if i = 0, j > 

ai(p,C)bo(C,p) if i > 1, J = 
ai(p,C)bj(C,p) 



(2^ 



Pi - c 



if « > 1, J > 1 
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Then for c > and y >0 we have 



E 



v(C,c) X AxY{p,y)dy. 



Proof. We start with the formula from Lemma 1 in [1] 



E 



-qT^-z(X +-c) 



E 



-(^c(,)-c)j(Xe(,) > C) 



E 



Using formula (19) we find that 



E 



g-2(Xe(9) 



i>l 



Next, using (18) and (13) after some algebraic manipulations we find that 



[z + 0)E 



nr 



^ n>l Pn 



^ Pj - Pj + ^ 



Combining the above three equations we find 



E 



-qT+-z{X +-c) 



5^a.(p,C)bo(C,p)Oe-^'^ + EE 



a.(p,C)bj(C,p) ^^-C.c Pj 



i>i i> 
which allows for a straightforward inversion in z, thereby completing the proof. 



(29) 



(30) 



□ 



5 Entrance/exit problems for a finite interval 

The two sided exit problem is a long standing problem of interest from the theory of Levy processes and 
there are very few cases where explicit identities have been be extracted for the (discounted) overshoot 
distribution on either side of the interval. Classically the only examples with discounting which have 
been analytically tractable are those of jump diffusion processes with (hyper)-exponentially distributed 
jumps (see for example [26, 44]) and, up to knowing the so-called scale function, spectrally negative 
processes, [34]. Otherwise, the only other examples, without discounting, are those of stable processes, 
[41], and Lamperti-stable processes, [13]. 
Recall that 

r+ = inf{t > : Xt > a}, = mi{t > : < 0}. 

For / e £oo(IR), the space of positive, measurable and uniformly bounded functions on M, and x G M, 
define the following operators 



igf){x) = I(0<x<a)E. [e-^-"V(X,+ )l(r+<ro-) 
(^/)(x) = 1(0 < X < a)E Je-^-o7(x )I (ro- < r^) 



(31) 
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and for any subset / G M define 



iVif)ix) 



I{x G /)E^ 



-qro 



(32) 



These are bounded operators from £00 (I^) into £00 (I^)- Note also that if Xq < then Tq = and 



similarly if Xq > a then r+ = 0. For example, when x > a, 



-IT I 



'fix. 



f{x) and when x < 0, 



E, 



Following similar reasoning to Rogozin [41] (see also [25]), an application of the Markov property 
tells us that for all x G [0, a] we have 



E, 



E, 



+E, 



e-''-V(X,+ )I (r+<ro-) 
e-'^oj (r- < r+) Ex _ 



Thus we have the following operator identities: 



(33) 



'P[0,a] = G + GV^^oofi], 
"^[O.a] =G + G'P[a,oo)- 

Consider a space H of bounded linear operators on £00 (1^) with the norm 

11^11= sup WGfWoo- 

ll/l|cc<l 

In the next theorem, which holds for general Levy processes, we give a series representation for G in 
terms of Vio^a], "^(-00,0]) 'P[a,oo), "^(-00,0] and a similar one for G- In the case of stable processes and general 
Levy processes respectively, Rogozin [41] and Kadankov and Kadankova [25] have shown similar series 
representations. The presentation we give above only differs in that it takes the form of linear operators. 
This turns out to be more convenient for the particular application to the M-class of processes. 

Theorem 4. Assume g > 0. ^ pair of operators {G,G) given by (31) is the unique solution in to the 
system of equations (33). Moreover, this solution can also he represented in series form 



a, 00 

)V( 

—00, 0] + • 



(34) 



and 



G = 'P[0,a] — 'P[o,a]'P[a,oo) + 'P[o,a]'P(-oo,0]'P[a,oo) — 'P[0,a]'P[a,oo)'P(-oo,0]'P[a,oo) + 

where, in both cases, convergence is exponential in the following sense. There exist t] G (0, 1) and a 
constant C > such that if Sn is the sum of the first n terms in the series describing G then \ \G — Sn\ \ < 
Cr]'"'^^, with a similar statement holding for G- 

Proof. It suffices to prove the result for G, the proof for G follows by duality. For uniqueness, we note 
first that the norm of operators ^(-00,0] is stricly less than one: 



IP, 



(-00,0] I 



sup 



supl(a; < 0)E^. 



^ fix. 



< snpl{x < 0)E, 



En 



< 1 
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and similarly we find HPfa.oo)!! < 1- Uniqueness also follows from the latter fact. Indeed if we assume 
that we have another pair of solutions {G, G), then using the system of equations (33) we find that the 
difference Q — G must satisfy the equation 



Thus we find 



G - G = {Q - G)V[a,oo)P{^-oofl]- 



\g-G\\ < \\G-G\\ X ||P[,,oo)|| X ||P(_oo,0]| 



and since ||P[a,oo)|| V | |'P(-oo,o] 1 1 < 1 we conclude that ||^ — G|| =0 and G = Q. To establish the series 
representation (34) we use (33) to find 

Q = T-'[Q,a\ " '^[0,a]'?'{-oo,0] + QV{a,oo)V {-cyo<d\ 

By iterating this equation we obtain series representation (34), which converges at an exponential rate 
as described in the statement of the theorem with 77 = HPfa.oo)!! V ||P(-oo,o]|| ^ (0, 1). □ 



The next theorem converts the above general setting to the specific setting of M-processes. 

Theorem 5. (First exit from a finite interval) Let a > and define a matrix B = B(p,C,a) 
{Kj}i,j>o with 



hi 







if i = 0, j > 1 
if i > 0, j = 



(35) 



. Pi + 

and similarly B = B(p, (, a). There exist matrices Ci, C2 and Ci, C2 such that for x G (0, a) we have 



e-^^h {X^, e dy ; r+ < r,) 
e-^-o'i (x- e dy ; r, < 



v(C, a — xY X Ci + v(C, x)^ X C2 X v(p, y — a)dy 
v(C, xf X Ci + v(C, a-xfx C2I X v(p, -y)dy (36) 



These matrices satisfy the following system of linear equations 



Ci 
C2 



A - C2BA 
-CiBA 



Ci 
C2 



A - C2BA 
-CiBA 



(37) 



where A and A are defined by (28). The system of linear equations (37) can he solved iteratively with 
exponential convergence with respect to the matrix norm \ \ ■ W^q, where for any square matrix M, 



IMIloo = max"^ |M, 

III i>0 ^-^ 



1,3 \- 



(38) 



i>o 



Before proceeding to the proof of the above theorem, there is one technical lemma we must first 
address. 
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Lemma 3. Let H = B x A. Then ||H||oo < 1- 

Proof. Using definitions (28) and (35) of A and B we find that for i,j>l 

hofi 



bo(C,p)$^afc(p,C)ae-"^'= 

k>l 



k>l 



CkC^ 

pj - Ck 



hi,o = piho{C,p)^akip,C) 



Cke' 



fc>l 



hij = p,bj(C,p)^afc(p,C) 



k>l 



Ck + pi 

(pj -Cfc)(Cfc + Pi)' 



Our first goal is to prove tliat hij are positive for all i,j > (note that hi^Q > for all i > 0). Let 
us consider a random variable ^ = e(pj) + Xe(g), where as usual e{pj) is an exponentially distributed 
random variable with rate pj, independent of X and e(g). The Laplace transform of ^ is given by 



Pj TT ^ + ^ 



nr 



Pj+Z\^l + f' 

■' n>l Sn 



and using the partial fraction decomposition (12), after some algebraic manipulations we find that ^ has 
a probability density function 

-P(e<x)=p, J]a,(p,C)^^. 

Pj - Cfc 



k>l 



This proves that Hqj are all positive. To prove that hi are positive one can use identity 



E[e-^»(«-)l(e>a)] =p, J]a,(p,C) 



Cke 



-a.(k 



k>l 



(pJ -Ck)iCk + Pi)' 



Next we need to prove that | |H| |oo < 1- Define a row vector hj as the i-th row of the matrix H. Using 
Theorem 3 we check that 



E 



e I(X^+ - a G dy)\Xo = -rj, = h^ x v{p,y)dy 



where r/o = and r/j is independent and exponentially distributed with parameter pj. Thus we have 



i>o 



< E 



Xn = 



< 1, 



and using the already established fact that /i,,- are positive we conclude that ||H||oo < 1- 



□ 



Proof of Theorem 5. The operator Vi is an integral operator, with kernel 

l{x G /)E^ e"'^^°^I(X^+ G dy) =I{x e I,y > a)v{(, a - x)'^ x A x v(p, y - a)dy, 



17 



(see Theorem 3). We also have a similar formula for Vi. Using the infinite series representation (34) we 
find that there exist some matrices of coefficients Ci, C2 and Ci, C2 so that integral kernels of operators 
Q and Q can be represented in the form (36). Matrix equations (37) follow from operator identities (33) 
and the fact that 

j v(p, z) X v(C, a + z^dz = B. 


Equations (37) can be solved iteratively as follows: 

Ci"+^) = A + CS"^BABA, Cf) = 

c("+i) = a + cS")baba, cS°) = o, 

and the convergence C["^ — t- Ci and CS"^ — )■ Ci is exponential because of Lemma 3. Once we find Ci 
and Ci we find C2 = -CiBA and C2 = -CiBA. □ 



Remark 3. Assume that within the M-class the positive (negative) jumps come from a mixture of 
N < 00 {N < 00) exponential distributions. Then matrices B and B have size (A^ + 1) x (A^ + 1) and 
(A^ + 1) X (iV + 1), while matrices Ci and Ci have sizes (A^ + 1) x (iV + 1) and (iV + 1) x (iV + 1) and 
can be found by the relations 



Ci = a(i- baba) \ Ci = a(i- baba) 



The idea of writing down a pair of simultaneous equations (31) also works when considering the 
problem of first entrance into an interval. This problem has been considered earlier by [41] for the 
setting of stable processes and [26] for two sided Levy processes whose upwards jumps are exponentially 
distributed. 

Theorem 6 (First entrance into a finite interval). Let a > and define r = inf{t > : Xt G (0,a)}. 
Define a matrix M = B(^, p, a)^ and similarly M = p, a)'^ (see (35)). There exist matrices Ni, N2 
and Ni, N2 such that for y G (0, a) we have 

fv(C,-a;)^ X [Ni X v(p,y) + N2 X v(p,a-y)] X dy, x<0 

[e-'^^I {Xr G dy)] = < . ^ I- - - - 1 (39) 

\v{C,x — ayx Ni X v(/), a — y) + N2 X v(p, y) x dy, x > a. 

These matrices satisfy the following system of linear equations 

= A + AMN2 TNi =A + AMN2 , , 

= AMNi [N2 =AMNi, ^ ^ 

where A and A are defined by (28). This system of linear equations (40) can be solved iteratively with 
exponential convergence with respect to the matrix norm \ \ ■ ||oo 

The proof of this Theorem is very similar to the proof of Theorem 5, and we leave the details to the 
reader. 
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6 Ladder processes 



In this Section we derive several results related to the Laplace exponent z) of the bivariate ladder 
process {L,H). We are interested in numerical evaluation of this object since it is the key to many 
important fluctuation identities, such as the quintuple law at the first passage which was introduced in 
[16]. 

Everywhere in this section we will denote the characteristic exponent of the Levy process X by 
'if{z) = — lnE[exp(i2;Xi)]. Note, that the Laplace exponent ip^z) can be expressed in terms of the 
characteristic exponent as ip{z) = — \E'(— iz), this fact easily follows from (1). 

The following Theorem, which holds for general Levy processes, is a generalization of the expression 
for the Wiener-Hopf factors which can be found in Lemma 4.2 in [37]. 



Theorem 7. Assume that for some ei > and for all \z\ ^K. large enough we have 

1^(^)1 > c\z\'\ 

for some c> 0, and 



£2 



dz, 

exists for some (and then for all) e2 > . Then for Re(g) > and Re(z) > we have 

1 f f\n{q + '^{u)) ln(l + ^(u))~ 



(41) 



(42) 



K{q, z) = exp 



27ri J \ u — iz 



u 



du 



(43) 



Proof. We start with the following integral representation for K,{q,z) (see Theorem 6.16 in [34]) 



K(g, z) = exp 




-qt—zx 



) -P(Xt e dx)dt 



Assume that the process Xt has a non-zero Gaussian component (a 7^ 0), then F{Xt e dx) has a density 
which can be obtained as the inverse Fourier transform of the right hand side of (1). Using this fact and 
assuming that 7 > we obtain 




(e-*-^- - e-«*-"") -F{Xt G dx)dt 



271 




(e-t-7- 



-qt—zx 



) J j e-**(")-'""dMda;dt. 



(44) 



Applying Fubini Theorem and performing integration in x we find that the above integral is equal to 



27r 



dt /e-*(i+*H) 
t J V 7 + 



1 

2^ 



dt 



du 



z + iu 

-*{!+*(«)) 



7 + iw 



1 

-{ 

z + m 



-tiq+'fiu)) 



(45) 
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where in the last step we have used the fact that for 7 > and Re(z) > 



1 



1 



7 + iu z + iu 



du = 0, 



(which can be easily verified by residue calculus). Next, we apply Fubini Theorem and Frullani integral 
to the last integral in (45) and combining this result with (44) we obtain 




(e-*-7- _ e""*---) _p(Xt G dx)dt 



/ln(g + *(M)) \n{l + <if{u)) 



2m J \ u ~iz 



u — 17 



du. 



(46) 



Next we need to take the limit of the above identity as 7 — )■ 0+. For u small the integrand in the right 
hand side of (46) is bounded by 



C 



ln(l + \^{u) 



\u\ 



uniformly in 7 G [0, 7] for each 7 > 0, and for u large it can be bounded by 



ln(g + ^(u)) ln(l + ^(u)) 



u — \z 



M — 17 



ln(g + -^{u)) In (} + 



1(^-7). ■ M ■ \ 

[U — l'^)[U — IZ) M — 17 

0{\n{\m{u)\)u-^) + 0{{u'^{u)y^), 



again uniformly in 7 G [0,7]. Using the above two estimates and assumption (41) we conclude that the 
integrand in the right hand side of (46) is bounded by 



^ln(l_+^HD 



\u\ 



(1 + H)-", 



for some C > uniformly in 7 G [0,7], and the above function is integrable on M due to assumption 
(42). Thus we can apply the Dominated Convergence Theorem on the right hand side of (46) together 
with the Monotone Convergence Theorem on the left hand side whilst taking limits as 7 — )■ 0+. To 
finish the proof we only need to take the limit a — t- 0+, which can be justified in exactly the same way 
using the Dominated Convergence Theorem on the right hand side of (46) and weak convergence on 
the right hand side of (46). Note in particular that the characteristic exponent of X is continuous in 
the Gaussian coefficient and hence, by the Continuity Theorem for Fourier transforms, the law of X is 
weakly continuous in the Gaussian coefficient. □ 



Assumption (42) is a very mild one. It is satisfied by all processes except those which have unusually 
heavy tails of the Levy measure. It is satisfied by all processes in M-class, and more generally, by all 
processes for which there exists an e > such that n(]R \ (— x,x)) = 0{x~'^) as x — > +00 (for example, 
all stable processes have this property). While condition (41) is more restrictive, one can see that it only 
excludes compound Poisson processes and some processes of bounded variation which are equal to the 
sum of its jumps. One example of such a process is a pure jump Variance Gamma process with no linear 
drift Xt = aWrt + fiTt, which satisfies |^E'(-2)| ~ cln(|2;|) as z — )■ 00. One can see that condition (41) 
is satisfied for the processes of bounded variation with non-zero drift, processes with non-zero Gaussian 
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component and all processes with the Levy measure satisfying n(]R \ {—x,x)) > as x — )■ O"*" for 
some e > 0. In particular this condition is satisfied for all processes in the 6-, hypergeometric or 
hyperexponential family excluding compound Poisson processes. 

Theorem 7 allows us to derive an expression for the Laplace exponent of the bivariate ladder process 
{L,H). 

Theorem 8. Assume that X belongs to /?-, 9-, hyper- exponential or hypergeometric family of Levy 
processes and that X is not a compound Poisson process. Then 

K(g,z)=[0,Mi.)]-^n77TT- ^4^) 

n>l ^'^^^^ 

Proof. If X belongs to /?-, 9- or hypergeometric family of Levy processes, we can use asymptotics for Cn as 
n ^ +CX) (see [29], [30] and [31]) to find that Cn(g)/Cn(l) = l + 0{n-^-') and Cn(g)/Cn(l) = l + 0(ra-^-^), 
which implies that both infinite products 

TT Cn{q) TT Cn(g) 

converge. Next, using the same technique as in the proof of Lemma 6 in [29], one can prove that 

Cn(g) TT ^ Cn(q) TT Cn(l) TT Cn(l) 
n>l ^ Cn(l) n>l ^ ^ C„(l) n>l ^''^^> n>l <^"WJ 

as z — )■ cxD, z G M. Using the Wiener-Hopf factorization q/{q + '^{z)) = (j)'^{z)(f)'^{z) and (18) we obtain 

q + ^{z) ^ Tj ^_^Asirj __Usl 

l-\-\lf(z) 111 11 1 -L ' 

^^^^^^ n>l^ Cn(l) n>l ^ + C„(l) 

Since X is not a compound Poisson process, we have |^E'(-2)| — )■ oo as z — )■ oo, ^ G M, thus we finally 
conclude that 

n>l ''"^^^ n>l Cn(lj 

The next step is to use the Wiener-Hopf factorization q/{q + "^{z)) = (f)g{z)(j)~{z), (18) and the above 
identity to rewrite q + "^{z) as 

n>l Pn n>l Pn 

Now we can use the above factorization and the following integral identity (which can be proved by 
shifting the contour of integration in the complex plane) 

1 f (\n{\{a-u)) \Yi{\{h-u))\ _ Jin (^(a - i^)) , if Im(a) < 0, Im(6) < 

(J. Uj 



2711 J \ u - iz u J 10, if Im(a) > 0, Im(6) > 0, 
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to deduce that for Iie{z) > 



1 
27ri 



ln(g + ^(M)) ln(l + ^(M))' 



u — IZ 



u 



In 1 + 



which is equivalent to (47). 



□ 



Corollary 4. Let z/(ds, dx) he the Levy measure of the ascending ladder process {L^ ^,Ht). Then for 
X > we have 



/e-«V(dM.) 



HCpV X v(p,x)dx. 



(4J 



Proof. Formula (48) is a corollary of (47), (17) and the fact that z/(ds,da;) is related to K,{q,z) through 
the formula 

oo oo 

K{q,z) = K{q,0) + az + J (l - e"^^) j e-'^'u{ds,dx). 



□ 



7 More Fluctuation identities 

We offer some more fluctuation identities. Although they are slightly more complex, they are still equally 
straightforward for the purpose of numerical work. 

We assume throughout this section that X is regular for both (0, oo) and (— oo, 0). Equivalently, we 
assume that ao(p, C) = ao(p, C) = 0. This is the case if, for example, X has paths of unbounded variation. 
It will be clear from the proofs of the results given below how this assumption may be removed. 

For a > and for y < a we define resolvent for X killed on leaving (— oo, a] as 

POO 

/?(")(«, dy) := / e-^*P(Xt e dy;t < r+)dt. 
Jo 

Theorem 9. Define a matrix D = {rfj,j}i,j>o follows 







if z = or j = 



^a,(p,C)-^a,(p,C) if^>l, J>1. 

Q + Cj 



Then if y < a we have 



qR^'i\a, dy) = v(C, V y)^ x D x v(C, V {-y))) - v(C, a)^ x D x v(C, a - y) 



dy. 
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Proof. From the proof of Theorem 20 on pl76 of [5], it can be seen that 



q j fiy)R^'^\a,dy) = y P(Xe(,) G d^) j P(-X,(^) G dM)/(-M + z) 

(-oo,a] [0,a] [0,00) 

P(Xe(,)Gd^) J /(y)P(X,(^)G-z + dy) 

[0,a] —(00,2] 

fiy) J P(Xe(,)Gdz)P(X,(^)G-^ + dy). 

(-oo,a] [OVj/,a] 

Thus we obtain an alternative representation of the Spitzer-Bertoin identity 

qR^'^\a, dy) = j P(Xe(,) G dz)P(Xe(,) e + d?/). 

[OVy,a] 

To finish the proof we have to use the formulas in (19) and perform the integration in the above expression 
(noting that some terms are lost on account of the fact that we have assumed ao(p, C) = ^o{p, C) = 0). □ 



Next define 

POO 

Q^'^\a,x,dy)= / e-^*P,(Xi Gdy,t<r+Aro-)dt 
Jo 

Theorem 10. Assume q > and y G [0,a], then 
qe^''\a,x,dy) 

v(C, OVy-xfxBx v(C, V (x - y)) - v(C, a - xf x B x v(C, a - y)] dy 
v(C, xf X Ci + v(C, a-xfx C2] [B{y) x D x v(C, 0) - B(a) x D x v(C, a - y)] dy. 

where matrix B(?/) = B(p,(,y) is defined in (35), D is defined in Theorem 9 while matrices Ci and C2 
come from Theorem 5. Also, v((^, 0) is intepreted as [0, Ci, C2) •••]'^- 



Proof. Define 



g^'^\a, X, dy) = e-'^o i(_x _ g dy ; r, < r+ 



(49) 



Note moreover that for / supported in [0, a] 

1, 



[0,a] 



f{y)Q'^''\a,x,dy) = -E, [/(X,(,))I(e(g) < r, A r+)] 



-E. [/(Xe(,))I(e(g) < r+)] - -E. [fiX^MMro < e{q) < r+)] 
r fiy)R^'\a -x,dy~x)- -E,. [/(Xe(,))I(ro- < e(g) < r+)] 

Jo q 
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The second expectation in the above expression can be rewritten as 



e-''*E, [f{Xt)I{T, < t < r+)] dt 



POD 

/ e-'^'fiX,)Iit<T:)dt 



g^'\a,x,dz) J R^'\a + z,z + dy)f{y). 



K+ 



So in conclusion we have for y E [0, a] 



e(«) (a, X, dy) = (a - X, dy - x) - / ^^^^ (a, x, dz)i?(«) {a + z,z + dy) 



and to end the proof one should use results of Theorems 5 and 9 and compute the above integral. The 
details are left to the reader. □ 



Remark 4. The previous result also allows us to write down the discounted joint overshoot, undershoot 
distribution for the two-sided exit problem: 

g^'^\a,x,dy,dz) = E, [e-'^h{T^ < ; G dy ; X^+_ G dz) 

(Again this relates to the so-called Gerber-Shiu measure for classical risk theory). Indeed, for ?/ > a and 
z G [0, a], using the compensation formula we have 



g^''\a,x,dy,dz) 



E, 



E. 



Y^e-'^Xt- < a, > 0, Xt^ G dz)I{Xt G dy) 



j e-«*I(Xt_ < a, X^_ > 0, G dz)dt 



U{dy - z) 



= e^''\a,x,dz)n{dy- z). 
Remark 5. Define the reflected process Y := X — 2L and its resolvent when killed on exiting [0, a] 

poo 

M'^\a, X, dy) = / e-'^'F,{Yt E dy, t < (T,)dt, 



where cTq = inf{s > : > a}. According to [4] one may write for all g > 0, 

A(^) (a, X, dy) = e('^) (a, x, dy) + g^^^ (a, x) ■ lim f . 

24-0 1 — g^'^>{a, z) 
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where g^'^\a,x) = Kx[e '^^o ^{Ta > Tq)]- Moreover, note that we may continue the computations in a 
similar way to before with the help of the compensation formula, 

Eo [e~'''^I{Y^^ G dy ; y;„_ G dz)] = M'^\a,0,dz)U{dy - z), 

and hence 

[e-'?'^"I(F.„ G dy ; G dz)] = g^''\a, x, dy, dz) + g^'^\a, x)M'^\a, 0, dz)U{dy - z). 



8 Numerical Results 



For all our numerical examples we will use a process X from the /3-family (see the introduction or [29]) 
having parameters 

(a, /i, ai, /3i, Ai, Ci, aa, ^2, A2, Ca) = (a, yU, 1, 1.5, 1.5, 1, 1, 1.5, 1.5, 1). 

Here /i = E[Xi] and a is the Gaussian coefficient, the other parameters define the density of a Levy 
measure, which has exponentially decaying tails and 0(|a;|~^/^) singularity at x = 0, thus this process 
has jumps of infinite activity but finite varation. We define the following four parameter sets 



Set 1: cr = 0.5, /i = 1 
Set 3: cr = 0, = 1 



Set 2: a = 0.5, /i = -1 
Set 4: (J = 0,/i = -1. 



As the first illustration of the efficiency of our algorithms we will compute the following three quan- 
tities related to the double exit problem: 

(i) density of the overshoot one the event that the process exists at the upper boundary 

d 



fi{x,y) 



dy 



E, 



(ii) probability of exiting from the interval [0, 1] at the upper boundary 



f2{x) 



E, 



(iii) probability of exiting the interval [0, 1] by creeping across the upper boundary 

Ux) = E, [e-'^-h [x^^ = 1 ; r+ < r^) . 

In order to compute these expressions we use methods described in Theorem 5. We truncate all the 
matrices Cj, Cj so that they have size 200 x 100 (this corresponds to truncating coefficients aj(p, Q and 
aj(/3, C) at 2 = 200 and coefficients bj(C,p) and bj(C,p) at j = 100). In order to compute coefficients 
a.i{p,C), ai(p, C), hj{(,p) and hj{(,p) we truncate infinite products in (14) and (15) at = 400, thus 
all the computations depend on precomputing {CnXn} for n = 1,2, ..,400. All the code was written in 
Fortran and the computations were performed on a standard laptop (Intel Core 2 Duo 2.5 GHz processor 
and 3 GB of RAM). 
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We present the results for g = 1 in Figures 1 and 2. Computations required to produce graphs for 
each parameter set took around 0.15 seconds. The numerical results clearly show the effects that we 
would expect to see. In Figures lb, Id and 2b we see a positive probability of creeping, which is expected 
since the process X has a Gaussian component in the first two cases and a bounded variation and positive 
drift in the third case. Parameter Set 4 corresponds to a process with bounded variation and negative 
drift, thus we do not have any upward creeping, and this is exactly what we obtain in Figure 2d. Also, 
figures la, Ic and 2a show that fi{x,y) — as x — )• 1~, which confirms our expectation, as in this case 
the upper half line is regular and as x — i- 1~ the process will cross the barrier at 1 by creeping, not by 
jumping over it. This is different from figure 2c, where, because of bounded variation and negative drift, 
the upper half line is irregular and the only way to cross the barrier at 1 is by jumping over it. Next, 
Figures lb and Id show that fi{x) — >• as x — t- 0+ and fi{x) — )■ 1 as x — 1~, which again agrees with 
the theory, as for parameter Sets 1 and 2 the process X has a Gaussian component, therefore is regular 
for {—oo, 0) and (0, oo). This is not so for parameter Sets 3 and 4, since now the process X has bounded 
variation and drift, and depending on the sign of the drift, is regular for either (0, oo) or (— oo, 0), and 
this is what we observe on figures 2b and 2d. 

For the next example, we compute the density u{s,x) of the bivariate renewal measure W(ds,dx) 
defined by 



where {L,H) is the ascending ladder process, see [16], [34] and [35]. This measure is a very important 
object, as it gives us full knowledge of the quintuple law at the first passage, see [16]. Using formulas 
6.18 and 7.10 in [34] we see that W(ds,dx) satisfies 



for any go > 0. It turns out that it is quite easy to compute the above integral numerically using 
technique discussed in [29] coupled with a Filon-type method (see [23]). Producing each graph on Figure 
3 takes around 1.2 seconds. In order to compute u{s,x) we truncated aj(p,C) at i = 100 and used 200 
roots (k to compute these coefficients using formulas (14). We chose go = 0.25, truncated integal in g 
at |Im(g)| < 10^. The numerical results are presented in Figure 3. Note that both the ascending ladder 
height and time processes behind Figure 3 (c) have a linear drift where as in Figure 3 (d) they are both 
driftless compound Poisson processes. In the former case this explains the strong concentration of mass 
around a linear trend, and in the latter case there exists an atom at x = s = 0, which is not visible on 
the graph since we are only plotting the absolutely continuous part. 

As we have mentioned in the introduction, all expressions related to fluctuation identities presented 
in this paper have the following property: (i) they are computed explicitly in terms of roots/poles of 
g + "^{iz) and possibly some linear algebra operations, (ii) all of them have the law of the space variables 
(for example, overshoot or location of the last maximum) in closed form and (iii) they involve Laplace 




R+ 




R+ 



therefore using Theorem 8 and (19) we find that the density of U{ds, dx) can be computed as 
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transform of the first passage time r+ or . The third condition imphes that if we want to compute 
joint distribution of both space and time functionals of the process (for example, joint density of the first 
passage time and the overshoot) we would have to perform a Fourier transform in the g- variable, and it 
has to be done numerically. See [29] and [31] for examples of application of this technique. It turns out, 
that using exactly the same method we can also be used to obtain similar reults for Levy processes with 
stochastic volatility, which are very popular models in Mathematical Finance, see for example [12]. We 
will briefly present this technique for the case of the two-sided exit problem considered above. Let Tg 
be an increasing continuous process satisfying Tq = 0. We require that E[exp(gTs)] is known in closed 
form, a classical example is the integral of the Cox-IngersoU-Ross diffusion process, however one can 
also choose several other processes, see [22]. Define a time-changed process Zg = Xt^, s > 0, where we 
assume that Tg is independent of Xt. As before, define Sq {s^} to be the first passage time of process 
Zg below {above a}. Since T, is continuous we have = r+, thus we obtain for any positive go 

[I edy; st <u; st < s^)] = [l (X^+ e dy ; r+ < r„ ; r+ < r,)] 
= Je^I edy;r^<t; r+ < r,)] P(T„ ^ dt) = ^ J g^'^\a, z, dy)E [e"^-] q-'dq, 

where we have used the fact that 

E, [I (X,+ e dy ; r+ < t ; r+ < r^)] = J 9^'^\a, z, dy)e'^'q-'dq, 

50+iM 

which follows from (49) by inverting the Laplace transform in q. Thus we see that if we choose the 
time change process {Tg : t > 0} for which the Laplace transform E[exp(gTs)] is known in closed form, 
then computing quantities for the time-changed process Z is essentially identical to computing the same 
quantities for the process X itself. 
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(a) Set 1: h{x,y) 



(b) Set 1: f2{x) (solid) and 73(2;) (circles) 




(c) Set 2: fi{x,y) (d) Set 2: f2{x) (solid) and /3(a;) (circles) 



Figure 1: Unbounded variation case (cr = 0.5): computing the density of the overshoot /i(x, y) (x G (0, 1), 
y G (0, 0.5)), probabihty of first exit f2{x) and probabihty of creeping fsi^x) for parameter Set 1 (positive 
drift /i = 1) and Set 2 (negative drift /i = — 1). 
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(a) Set 3: fi{x,y) (b) Set 3: f2{x) (solid) and 73(2;) (circles) 




(c) Set 4: fi{x,y) (d) Set 4: f2{x) (solid) and /3(a;) (circles) 



Figure 2: Bounded variation case (cr = 0): computing the density of the overshoot fi{x,y) (x G (0, 1), 
y G (0, 0.5)), probabihty of first exit f2{x) and probabihty of creeping fsi^x) for parameter Set 3 (positive 
drift /i = 1) and Set 4 (negative drift /i = — 1). 
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(c) Set 3 (d) Set 4 

Figure 3: Computing the density of the renewal measure U{ds, dx), s G (0, 0.25) and x G (0, 0. 
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